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Since  current  empirical  models  specifying  low-latitude  electron  density  profiles  severely  underesti¬ 
mate  the  daytime  plasma  density  scale-height  and  total  electron  content  (TEC)  values,  a  semi-empirical 
low-latitude  ionospheric  model  (SLIM)  has  been  developed  which  is  not  only  computationally  fast,  but 
also  more  realistic.  Electron  density  profiles  (180  1800  km)  are  theoretically  calculated  as  a  function  of 
latitude  (every  2  between  24  N  and  24  S  dip  latitude)  and  local  time  levery  half  hour,  over  24  hours 
LT)  by  solving  the  time-dependent  plasma  continuity  equation  Using  simple  exponential  functions, 
sets  of  coefficients  are  then  generated  which  reproduce  these  individual  profiles  The  coefficients  them¬ 
selves  are  easily  stored,  quickly  retrieved  and  form  the  basis  for  a  fast,  portable,  semi-empirical  com¬ 
puter  code.  This  paper  describes  briefly  the  input  parameters  used  to  theoretically  calculate  the  profiles 
and  the  procedures  used  to  generate  the  coefficients.  The  SLIM  profiles  are  compared  with  the  Chiu 
and  Bent  empirical  models  for  equinox,  solar  maximum  conditions,  while  calculated  6100  A  airglow 
intensities  and  TEC  values  are  compared  with  available  observations  The  SLIM  profiles,  their  coef¬ 
ficients.  TEC  and  6300  A  airglow  intensities  are  available  in  tabular  and  computer  formats. 


I.  INTRODUCTION 

1 . 1 .  Background 

Current  empirical  models  of  Ihe  low-latitude  iono¬ 
spheric  F  region  [Llewellyn  and  Bent,  1973;  Chiu. 
1975;  Rawer.  1981]  can  severely  underestimate  the 
daytime  plasma  density  scale  height  and  total  elec¬ 
tron  content  (TEC)  when  compared  with  actual  ob¬ 
servations.  As  illustrated  in  Figure  1,  the  Internation¬ 
al  Reference  Ionosphere  (IRI)  described  by  Rawer 
[1981]  yields  daytime  TEC  values  which  are  50% 
lower  than  observed  TEC  values  at  Manila  (dip 
14.5  N)  in  January  1982  [ McNamara ,  1984].  Some 
improvement  occurs  when  the  Bent  topside  model  is 
incorporated  with  the  IRI  bottomside  models  but  the 
predicted  values  are  still  lower  than  observed.  Sub¬ 
stantial  improvement  is  achieved  when  theoretically 
calculated  profiles  [Anderson.  1973]  are  used  to  pre¬ 
dict  TEC  values.  The  reason  is  that  vertical  plasma 
transport  by  upward  E  x  B  drift  produces  both  top¬ 


side  and  bottomside  profiles  which  are  much  broader 
(thicker)  than  Chapman-like  profiles. 

To  theoretically  calculate  electron  density  profiles 
whenever  they  are  needed,  however,  is  prohibitively 
time-consuming  on  even  the  fastest  computer.  An 
alternate  solution  to  this  problem  is  to  theoretically 
calculate  electron  density  profiles  as  a  function  of 
latitude  and  local  time  and  then  generate  coefficients 
which  produce  these  profiles.  The  coefficients  them¬ 
selves  are  easily  stored,  quickly  retrieved  and  form 
the  basis  for  a  fast,  portable  computer  code  which 
will  produce  realistic  low-latitude  F  region  electron 
density  profiles.  The  results  of  this  technique  are 
called  the  semi-cmpirical  low-latitude  ionospheric 
model  (SLIM). 

Using  the  techniques  described  by  Anderson 
[1973]  electron  densities  as  a  function  of  altitude, 
latitude  and  local  time  are  calculated  by  solving  the 
time-dependent  ion  (O*)  continuity  equation  nu¬ 
merically.  The  effects  of  production  by  photoioniza¬ 
tion,  loss  through  charge  exchange  with  N2  and  02 
and  transport  by  diffusion,  neutral  winds  and  vertical 
E  x  B  drift  are  included.  Because  E  x  B  drifts  play 
such  an  important  role  in  determining  the  electron 
density  distribution  [Hanson  and  Moffett.  1966;  An- 
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t  tg.  1  A  comparison  between  observed.  predicted,  and  theo- 
reticall)  calculated  total  electron  content  tTI  Cl  values  at  Manila 
as  a  function  of  local  time  for  January  19X2.  R  is  the  sunspot 
number  and  DIP  refers  to  the  angle  the  geomagnetic  held  line 
makes  with  the  local  hori/on. 

derson.  1973]  and  have  a  seasonal  and  solar  cycle 
dependence  [ Fejer  et  al..  1979],  the  calculations  are 
carried  out  for  equinox  and  solstice  conditions  for 
both  solar  cycle  maximum  and  minimum  periods. 

Once  electron  density  profiles  arc  obtained,  they 
are  normalized  to  the  peak  electron  density, 
and  appropriate  coefficients  which  reproduce  the 
normalized  profiles  (180  1800  km)  every  half  hour 
local  time  and  every  2  dip  latitude  are  calculated. 
To  regenerate  density  profiles  simply  requires  speci¬ 
fying  Nmat  internally  (that  is,  reuse  of  the  theoretical 
.Vm„  values)  or  by  externally  using  an  empirical 
model  such  as  CCIR  [1978]  or  observed  faF2  values, 
and  then  choosing  the  calculated  coefficients  appli¬ 
cable  to  the  specified  latitude  and  local  time  location. 
In  the  following  sections,  we  describe  briefly  the 
various  input  parameters  which  were  used  in  calcu¬ 
lating  theoretically  the  electron  density  profiles  from 
the  time-dependent  continuity  equation.  In  the  last 
sections  we  compare  (1)  the  current  Chiu  and  Bent 


empirical  profiles  with  SLIM  profiles  at  specific  lati¬ 
tudes  and  local  times  for  equinox,  solar  cycle  maxi¬ 
mum  conditions.  (2)  calculated  and  observed  6300-A 
airglow  intensities,  and  (3)  calculated  and  observed 
TEC  values. 

1.2.  Model  calculations 
The  ion  (O  * )  continuity  equation  is  given  by 

r.V 

— -  +  V  .</V,V ,)  =  /»,-  /,,  (I) 

cl 

where  =  ,Vt.)  is  the  ion  density;  P, ,  the  ion  pro¬ 
duction  rate;  the  ion  loss  rate;  and  V,,  the  ion 
transport  velocity.  In  the  ionosphere,  plasma  is  trans¬ 
ported  along  geomagnetic  field  lines  by  diffusion  ami 
neutral  winds  and  perpendicular  to  field  lines  pri¬ 
marily  by  E  x  B  drift  [Kendall  and  Pickering.  1967]. 
Solving  (I)  at  low  latitudes  is  facilitated  by  trans¬ 
forming  the  independent  coordinates  r.  1 1  and  4>  to  a 
coordinate  system  parallel  and  perpendicular  to  B 
[see  Anderson.  1973],  Equation  (I)  can  be  written 

'‘•V, 

—  +  V,  •  V.v  -  P,  -  -  V  •  (.V|  V,  )  -  V  •  V,  (2) 

cl 

where  V,  is  given  by  E  x  B  B2  and  V,  includes  the 
effects  of  plasma  diffusion  and  neutral  wind.  The 
right-hand  side  of  the  equation  involves  terms  which 
are  second  order  in  the  coordinate  parallel  to  B.  The 
left-hand  side  of  (2)  is  the  time  rate  of  change  of 
plasma  density  in  a  frame  of  reference  which  drifts 


Eig.  2.  Ion-neutral  and  eicciron-neutral  temperature  ratios  as 
a  function  of  altitude  The  neutral  temperature  is  adopted  from 
the  MSIS  model  atmosphere  (see  text  for  details). 
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Eig.  3 b.  Vertical  E  x  B  plasma  drift  as  a  function  of  local  time 
during  solar  cycle  maximum  and  minimum  periods  for  June  sol¬ 
stice  conditions. 


with  the  E  x  B  drift  velocity.  Equation  (2)  is  solved 
numerically  to  give  ion  densities  as  a  function  of  alti¬ 
tude,  latitude  and  local  time  (see  Moffett  [1979]  for  a 
review  of  transformations  and  numerical  solutions). 

The  set  of  coefficients  for  the  ion  continuity  equa¬ 
tion  is  obtained  from  models  of  the  neutral  compo¬ 
sition,  ion  and  electron  temperatures,  and  pro¬ 
duction,  loss  and  diffusion  rates  as  well  as  E  x  B 
drift  and  neutral  wind  models.  Briefly,  the  models  are 
the  following: 

1.  The  MSIS  [ Hedin ,  1983]  neutral  atmospheric 
model  is  used  to  calculate  N2,  02,  and  O  densities 
and  the  neutral  temperature,  T„,  as  a  function  of  alti¬ 
tude,  latitude  and  local  time.  An  El 0.7  cm  flux  of  180 
units  is  chosen  to  represent  solar  cycle  maximum 
period,  and  70  units  for  solar  minimum  conditions. 
Average  geomagnetic  conditions  were  chosen  by  set¬ 
ting  Ap  =  15  in  the  MSIS  model. 

2.  Production  and  diffusion  rates  are  similar  to 
those  used  by  Anderson  [1973],  For  the  photoioniza¬ 
tion  rate  at  the  top  of  the  atmosphere,  P  r  ,  a  value  of 
5.5  x  10"’  s  1  represents  solar  cycle  maximum  con¬ 
ditions,  with  a  value  of  2.3  x  10  7  s'1,  for  solar 
minimum  conditions.  The  loss  rate  coefficients  are 
adopted  from  Torr  and  Torr  [1979], 

3.  Plasma  temperatures  were  chosen  by  adopting 
TJT„  and  7j/T„  ratios  as  a  function  of  local  time, 
altitude  and  season  from  the  IRI  model  and  then 
applying  those  ratios  to  the  T„  results  obtained  from 
MSIS.  Figure  2  illustrates  the  height  variation  of  the 
ratios  for  the  sunspot  maximum  equinox  case. 

4.  The  vertical  E  x  B  drift  velocity  observed  by 
Woodman  [1970]  is  essential  in  producing  and  main¬ 
taining  the  electron  density  distribution  observed 
near  the  magnetic  equator.  We  adopt  six  vertical 


VERTICAL  E«e  DRIFT 


f-  ig.  V  Vertical  E  x  B  plasma  drift  a«  a  function  of  local  time 
during  solar  cycle  maximum  and  minimum  periods  for  equinox 
conditions. 


drift  models  to  represent,  respectively,  three  seasons 
(equinox,  summer  and  winter  solstices)  at  solar  cycle 
maximum  and  solar  cycle  minimum  periods.  The  dif¬ 
ferent  drift  patterns  arc  based  primarily  on  the  Jica- 
marca  incoherent  scatter  radar  observations  reported 
by  Fejer  et  at.  [  1 979],  We  have  assumed  the  vertical 
drift  velocity  is  independent  of  altitude.  The  east-west 
component  of  E  x  B  drift  reported  by  Fejer  et  al. 
[1981]  is  neglected  in  the  calculations  because  its 
effect  on  electron  density  profiles  is  negligible  [An¬ 
derson,  1981].  Figures  3a,  3b  and  3r  display  the  diur¬ 
nal  variation  in  drift  for  the  equinox,  June  solstice 
and  December  solstice  cases,  for  solar  maximum  and 
minimum,  respectively. 

5.  The  geomagnetic  field  is  represented  by  a 
dipole  model  whose  axis  is  aligned  with  the  earth’s 
rotational  axis.  Although  we  incorporate  E  x  B  drift 
models  appropriate  to  the  American  sector  (where 
the  dip  equator  lies  12  south  of  the  geographic 
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Pig.  m  .  Vertical  t  x  B  plasma  drift  as  a  function  of  loc«l  time  <3q9 
during  solar  cycle  maximum  and  minimum  periods  for  December  ^ 
solstice  conditions. 
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LOCAL  TIME  (hrs) 


Fig.  4u.  Concour  plot  of  meridional  neutral  wind  as  a  function  of  latitude  and  local  lime  for  equinox  periods 

Positive  values  represent  north-to-south  winds. 


equator),  longitude  effects  are  not  addressed  at  this 
lime.  Further  development  of  the  model  will  include 
a  tilted  dipole  geomagnetic  field  with  declination 
specified  over  the  regions  of  interest. 

6.  A  simple  analytic  formulation  for  the  meridio¬ 


nal  neutral  wind  field  is  adopted  which  makes  use  of 
observations  and  theory.  It  is  constructed  on  the 
basis  of  observed  mean  (DC),  diurnal  (24-hour),  semi¬ 
diurnal  (12-hour),  and  terdiurnal  (8-hour)  harmonics 
at  Arecibo  (18  N).  Millstone  Hill  (42  N)  and  St. 


JUNE  SOLSTICE  MERIDIONAL  WINDS  (m/s) 


Fig  4ft.  Contour  plot  of  meridional  neutral  wind  as  a  function  of  latitude  and  local  time  for  June  solstice  periods. 

Positive  values  represent  north-to-south  winds 
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DECEMBER  SOLSTICE  MERIDIONAL  WINDS  (m/s) 


Kig.  4<  Contour  plot  of  meridional  neutral  wind  as  a  function  of  latitude  and  local  time  for  December  solstice 
periods.  Positive  values  represent  north-to-south  winds. 


Santin  (45  N).  with  latitude  variations  of  amplitude 
and  phase  obtained  from  theoretical  (numerical)  sim¬ 
ulations  [Forbes.  1982a.  /»].  The  general  latitudinal 
pattern  assumed  is  guided  by  the  tidal  calculations  of 
Forbes  [1982a.  h]  and  the  zonal  mean  velocities  pre¬ 
sented  by  Roble  et  al.  [1977],  The  analytic  ex¬ 
pressions,  valid  only  for  latitudes  |  (1|  <45  ,  are  given 
in  the  appendix  for  equinox.  December  solstice  and 
June  solstice  periods. 

For  solstice  conditions,  the  diurnal  component, 
F24.  is  shifted  10  in  latitude  in  the  direction  of  the 
solar  bulge  (consistent  with  Forbes).  The  Roble  el  al. 
[1977]  mean  meridional  winds.  V0.  during  solstice 
are  shifted  downwards  by  about  40  m/s  so  that  the 
transequatorial  amplitudes  are  of  the  order  50  m/s 
and  the  winter  hemisphere  reversal  occurs  near  Mill¬ 
stone  Hill. 

The  solstitial  modifications  to  the  mean  and  diur¬ 
nal  components  adequately  account  for  the  main  fea¬ 
tures  of  the  observed  seasonal-latitudinal  dependence 
of  the  meridional  wind  for  our  purposes.  Seasonal 
variations  in  the  semidiurnal  and  terdiurnal  compo¬ 
nents  are  neglected.  Figures  4a,  4b.  and  4c  give  con¬ 
tour  plots  of  the  overall  wind  model  spanning  the 
simulation  region  where  positive  values  represent 
north-to-south  winds.  Finally,  it  should  be  noted  that 
while  EUV  varies  markedly  over  a  solar  cycle,  the 
associated  variations  in  ion  drag  suppress  solar  cycle 


variations  in  the  tidal  winds,  although  neutral  tem¬ 
peratures  vary  considerably.  Therefore,  we  have  as¬ 
sumed  the  current  wind  model  parameterization  to 
be  independent  of  solar  activity. 

2  THh  SKMI-EMP1RICAL  APPROACH 

2.1.  The  SLI M  dam  base 

As  described  in  section  1.1,  the  theoretical  model  is 
used  to  produce  six  preliminary  data  bases.  Three  of 
these  runs  are  for  sunspot  minimum  (F10.7  =  70, 
Ap  —  15)  during  the  June  solstice,  December  solstice, 
and  equinox  seasons.  The  remaining  three  are  for 
sunspot  maximum  (F  10.7  =  1 80,  AP  =  15)  during 
the  same  seasons.  Each  record  of  data  consists  of 
profile  parameters  stored  for  a  fixed  geographic  (geo¬ 
magnetic)  point  and  fixed  time.  Each  data  base 
covers  25  latitudes  (24  S-24  N.  every  2°)  and  48 
times  (0000-2330  LT,  every  30  min).  A  subset  of  the 
SLIM  data  set,  using  a  more  coarse  resolution  for 
the  electron  density  profiles  (Ne(h)  every  4°  in  latitude 
and  every  hour  of  local  time)  was  presented  in  tabu¬ 
lar  format  by  Anderson  et  al.  [1985]. 

The  core  of  the  SLIM  approach  is  to  identify  a 
relatively  simple  Nf(h)  function  that  will  have  a  small 
number  of  coefficients  and  yet  provide  a  good  repre¬ 
sentation  of  the  profiles  calculated  from  the  theoreti¬ 
cal  model.  As  guiding  principles,  it  was  decided  that 
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ELECTRON  OENSlT  *  (n/cm'l 

Fig.  5.  A  comparison  between  SLIM-generated  electron  density  profiles  (open  circles)  and  the  empirical  models 
of  Chiu  (solid  squares)  and  Bent  (crosses)  for  equinox,  solar  cycle  maximum  conditions  at  1400  LT.  Profiles  at  the 
magnetic  equator  (top  portion)  and  at  16  S  dip  latitude  (bottom  portion). 

(1)  maximum  accuracy  was  required  near  the  F  days  of  ionospheric  physics  offer  p  simple  way  to 
region  peak  (in  both  density  and  height),  and  (2)  explain  the  vertical  structure  of  plasma  in  the  upper 
while  some  errors  in  Nt  could  be  accepted  well  below  atmosphere.  They  are  based  on  relatively  simple 
or  above  the  peak  density  height  the  overall  physics  and  thus  cannot  be  used  to  represent  the 

profile  (as  represented  by  the  total  electron  content  entire  Nf{h)  profiles  that  emerge  from  the  more  so- 
(TEC)  parameter)  had  to  reproduce  the  theoretical  phisticated  simulation  techniques  in  use  today.  As 
results  to  within  5%  or  less.  simple  mathematical  formulas,  however,  they  are  well 

The  Chapman  profiles  developed  in  the  earliest  suited  for  representing  portions  of  the  ionospheric 
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Fig  6.  A  comparison  between  SLIM-generated  electron  densily  profiles  lopen  circles)  and  Ihe  empirical  models 
of  Chiu  (solid  squares)  and  Bent  (crosses)  for  equinox,  solar  cycle  maximum  conditions  at  2000  LT.  Profiles  at  the 
magnetic  equator  (top  portion)  and  at  16  S  dip  latitude  (bottom  portion). 


profile.  Specifically,  while  a  single  Chapman  profile 
cannot  fit  the  entire  Ne(h)  profile,  separate  Chapman 
formulas  can  be  used  to  reproduce  the  topside  and 
bottomside  components  of  the  overall  profiles  with 
good  results.  It  was  decided,  therefore,  to  define  a 
“generalized"  or  “modified"  Chapman  formula  with 
coefficients  that  vary  for  each  topside  and  bottom- 
side  component  of  the  Ne(h)  profile.  The  function 


used  to  fit  the  theoretically  derived  profiles  was  thus 
chosen  to  be 

N,(h)  =  JV  exp  [dl  -  :  -  e'*)] 

(3) 

2  -  - 

A 

In  using  (3),  six  parameters  (or  coefficients)  are 
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ELECTRON  DENSITY  (el/ciYt1) 

Kig  7  A  comparison  between  SLIM-generated  electron  density  profiles  (open  circles)  and  the  empirical  models 
of  Chiu  (solid  squaresl  and  Bent  (crosses)  for  equinox,  solar  cycle  maximum  conditions  at  0500  LT.  Profiles  at  the 
magnetic  equator  (top  portion)  and  at  16  S  dip  latitude  (bottom  portion). 


needed  to  reproduce  the  A le(h)  profile:  Nmn. 
c‘Up .  -4up,  <  |0,  Alo.  The  advantages  of  this  method  are 
(1)  the  fit  achieves  100%  accuracy  at  and  hmtl 
since  they  are  themselves  input  coefficients:  (2)  the 
topside  and  bottomside  Chapman-type  functions 
guarantee  continuity  of  both  density  and  its  first  de¬ 
rivative  across  the  fit  boundary;  and  (3)  with  a  fit 
technique  weighted  by  density  values  (see  below),  the 
c  and  A  chosen  emphasize  maximum  accuracy  in 
regions  near  the  peak. 

Included  on  each  record  of  the  database  are  two 


other  numbers  available  for  use  and/or  inspection. 
These  are  (1)  the  total  electron  content  (TEC)  ob¬ 
tained  by  integrating  the  SLIM  profile  from  180  to 
1800  km;  and  (2)  6300-A  emission  calculated  from 
02+  recombination  chemistry  only  (and  therefore,  ap¬ 
propriate  for  postsunset  airglow). 

2.2.  Production  details 

The  report  by  Anderson  et  al.  [1985]  describes  the 
numerical  methods  employed  to  obtain  the  parame- 


|T 


ANDERSON  ET  AL  LOW-LATITUDE  IONOSPHERIC  MODEL 


300 


electron  density  <ei 1 » 


I  ig  8.  A  comparison  between  SI.IM-gencratcd  electron  density  profiles  (open  circles!  and  profiles  measured  by 
ihe  Jicamarcu  incoherent  scatter  radar  (solid  circles)  on  March  5.  196*).  at  1400  LT  (top  portion)  and  2000  LT 
(bottom  portion) 


ters  that  form  the  basis  for  SLIM.  Briefly.  \max  and 
arc  known  precisely  since  they  arc  taken  directly 
from  the  theoretically  derived  profiles:  the  coef¬ 
ficients  A  and  c  were  obtained  by  separate  least 
squares  fits  of  the  generalized  Chapman  formula  to 
the  electron  densities  above  and  below  the  peak. 

V  COMPARISON  WITH  EMPIRICAL  MODELS 

Figure  5  compares  the  semi-empirical  low-latitude 
ionospheric  model  (SLIM)  with  the  Chiu  and  Bent 
models  for  equinox,  solar  maximum  conditions  at 
1400  LT.  The  daytime  upward  E  x  B  drift  pictured 


in  Figure  3«  causes  the  broad  F  region  profile.  Re¬ 
spective  topside  scale  heights  at  700  km  are  540  km. 
150  km  and  1 15  km.  Integrations  of  the  electron  den¬ 
sity  profiles  with  altitude  give  total  electron  contents 
(TEC)  of  97.  48,  and  61  (  x  10' 2  el 'em2),  respectively. 
The  bottom  portion  of  Figure  5  shows  the  compari¬ 
son  at  16  S  dip  latitude  at  the  crest  of  the  equatorial 
anomaly.  Although  the  SLIM  topside  scale  height  is 
not  as  large  as  at  the  magnetic  equator,  it  is  still 
larger  than  those  given  by  Chiu  or  Bent.  In  addition, 
the  /Vmj<  value  of  3.5  x  I0h  cl/cm3  is  significantly 
larger  than  the  value  at  the  magnetic  equator  due  to 
the  combination  of  upward  and  poleward  E  x  B  drift 
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and  downward  diffusion.  Note  that  neither  the  Chiu 
or  Bent  profiles  shows  any  significant  change  with 
latitude. 

The  most  dramatic  differences  between  SLIM  and 
the  empirical  models  occur  just  after  sunset  around 
2000  LT  and  are  pictured  in  Figure  6.  At  the  mag¬ 
netic  equator  (top  portion),  the  postsunset  enhance¬ 
ment  in  upward  drift  lifts  the  F  layer  to  600  km 
compared  to  /imj(  values  of  450  km  and  350  km  for 
the  Bent  and  Chiu  models,  respectively.  Such  differ¬ 
ences  in  hmdx  would  cause  significant  differences  in 
the  estimated  6300-A  airglow  intensity  and  the  ion- 
neutral  drag  coefficient  in  the  neutral  momentum 
equation. 

At  16  S  dip  latitude  (bottom  portion  of  Figure  6) 
the  same  upward  drift  causes  an  enhancement  in 
/Vmal.  where  the  calculated  value  is  5  x  10*  el/cm1 
compared  with  2  x  10*  el  cm1  for  Chiu  and  8  x  10* 
el/cm3  for  the  Bent  model.  Respective  TEC  values 
(  x  I012  el/cm2)  are  139.  48.  and  19.  Finally,  Figure  7 
gives  a  comparison  between  the  SLIM,  Chiu  and 
Bent  at  0500  LT.  The  calculated  profiles  at  the  mag¬ 
netic  equator  and  16 'S  dip  latitude  are  much  more 
Chapman-like  and  in  reasonable  agreement  with  the 
empirical  models  in  both  hmmx  and  topside 

scale  height.  The  calculated  value  of  TEC  at  16°S  dip 
latitude  (bottom  portion)  is  5  x  1012  el/cm3  which 


compares  favorably  with  the  Beni  value  of  3.5  x  I012 
and  a  Chiu  value  of  9  x  I012  el  cm2. 

f  igure  8  compares  electron  density  profiles  mea¬ 
sured  on  March  5.  1969.  using  the  Jicamarca  inco¬ 
herent  scatter  radar  facility  with  calculated  density 
profiles  at  1400  LT  and  2000  LT.  The  comparison 
demonstrates  that  both  the  large  daytime  plasma 
scale  height  and  the  calcula'ed  values  of  /ima,  at  2000 
LT  are  in  reasonable  agreement  with  observed 
values.  Again,  it  should  be  emphasized  that  the  pro¬ 
files  are  sensitive  to  the  vertical  E  x  B  drfft  velocity 
whose  diurnal  variation  may  have  significant  day-to- 
day  changes. 

4  COMPARISON  WITH  OBSTRVHD  fttOO-A  AIRGLOW 
INTKNSITIFS 

Based  on  the  SLIM  density  profiles  and  the  MSIS 
neutral  atmosphere  for  the  equinox,  solar  cycle  maxi¬ 
mum  case,  we  calculated  6300-A  airglow  intensity  at 
night  using  the  following  formula: 

a  ~ 

rno("m  f _ *,[Q33[f  ] _ i  dh 

J i so  km  t  '4,ol  +  (k2[N2])  4-(t:  ,[e "  ]  +(k4[02])J 


(4) 
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Fig  10  Measured  63(X)  A  intensit)  contours  (Rayleighs)  as  a  function  of  latitude  and  local  time  under  equinox, 
solar  cycle  maximum  conditions  [after  ,4hreu  el  at..  1982]. 


where 

A,  -  2  x  10  "  (300  Tj'  ’  cnv' s  ' 

[Hanks  and  Km  karts.  1973]: 

k.  =  2.3  x  10  "  cm' s  1  [Torr  and  Ton.  1982]; 
ks  —  2  x  10  “  cm' s  1  [Sealon.  1955]: 

-  5.3  x  10  "  cm’ s  1  [ Ton  and  Ton.  1982]: 

•!„„  -  0.006K  I  s  1  [Ton  and  Torr.  1982]; 

=  0.00515  s’ 1  f Torr  and  Torr.  1982]; 
dh  -  height  element,  km. 

Figure  9  presents  6300-A  isointensity  contours  in 
Rayleighs  as  a  function  of  dip  latitude  arid  local  time. 
Near  the  magnetic  equator,  just  after  sunset,  the  in¬ 
tensity  drops  to  less  than  10  Rayleighs  due  to  the 
high  altitude  (over  350  km)  of  the  bottomside  of  the 
F  layer.  Peaks  in  airglow  intensity  occur  at  +16 
latitude  caused  by  the  crests  in  ionization  which 
occur  there  and  the  lower  altitude  of  the  F  layer.  As 
local  time  increases,  the  crests  in  airglow  intensity 
move  toward  the  magnetic  equator  due  to  downward 
E  x  B  drift  which  lowers  the  F  layer  as  it  transports 
plasma  back  toward  the  equator.  The  subsequent 
motion  of  the  airglow  crests  away  from  the  equator 
after  2300  LT  is  due  to  the  abatement  of  equator- 


ward  wind  after  2200  LT  which  allows  the  F  layer  to 
move  to  a  lower  altitude  thus  increasing  the  airglow 
intensity. 

F  igure  10  presents  iso-intensity  contours  of  6300-A 
airglow  based  on  data  collected  by  the  visible  air¬ 
glow  experiment  (VAH)  on  board  the  AE-E  satellite 
[Abri'ii  el  al„  1982].  These  observations  were 
averaged  over  longitude  from  September  1980  to  No¬ 
vember  1980  and  are  essentially  representative  of 
equinox,  solar  cycle  maximum  conditions.  Several  of 
the  features,  such  as  the  very  low  airglow  intensities 
near  the  magnetic  equator  after  sunset,  the  crests  in 
airglow  that  move  toward  the  magnetic  equator,  and 
the  subsequent  motion  away  from  the  equator  after 
2300  LT  are  reproduced  by  the  SLIM  calculated  air¬ 
glow  results.  The  absolute  intensity  levels  in  Ray¬ 
leighs  show  an  excellent  agreement  in  the  premid¬ 
night  sector:  the  postmidnight  enhancements  are  fac¬ 
tors  of  the  2-3  higher  in  SLIM  as  compared  to  the 
AF.  average  values  (an  efFect  due  to  neutral  wind  dif¬ 
ferences).  Since  6300-A  airglow  is  produced  near  the 
bottomside  of  the  F  layer  and  is  very  sensitive  to  the 
height  of  this  region,  the  basic  agreement  between 
calculated  and  observed  values  gives  us  confidence 
that  the  assumed  neutral  wind  model  and  the  E  x  B 
drift  model  are  fairly  realistic,  at  least  for  this  season 
and  solar  cycle.  A  paper  is  in  preparation  which 
examines  the  separate  and  combined  roles  played  by 
neutral  wind  and  E  x  B  drift  in  producing  the  ob- 
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SLIM  INDIAN  SECTOR 


LOCAL  TIME  (hrj| 

Fig.  II  Contours  of  calculated  TEC  values  <  x  10'-’  el  cm 2 1. 
for  equinox,  solar  cycle  minimum  conditions,  left:  and  contours  of 
observed  TEC  values  in  the  Indian  sector  [after  Klobuchar  el  al.. 
1977J,  right 

served  latitude  and  local  time  variations  in  6300-A 
intensity. 

5.  COMPARISON  BETWEEN  CALCULATED  AND 

OBSERVED  TOTAL  ELECTRON  CONTENT  VALUES 

Integration  of  an  electron  density  profile  with  alti¬ 
tude  gives  the  total  electron  content  (TEC).  Figures 
1 1  and  12  present  contour  plots  of  TEC  as  a  function 
of  latitude  and  local  time  calculated  from  the  equin¬ 
ox,  solar  cycle  minimum  and  maximum  SLIM  pro¬ 
files  (integrated  over  the  180-1800  km  altitude 
range).  Daytime  crests  in  TEC  values  occur  on  either 
side  of  the  equator  under  both  minimum  and  maxi¬ 
mum  solar  flux  conditions;  these  correspond  to  crests 
in  the  peak  electron  density,  Nmtx,  having  the  same 
latitudinal  distribution  (the  so-called  equatorial 
anomaly).  At  the  magnetic  equator,  however,  there  is 
only  one  daytime  maximum  value  in  TEC,  in  con¬ 
trast  to  the  two  maxima  characteristically  seen  in 
(one  >n  the  morning  and  one  in  the  afternoon). 
This  is  in  agreement  with  observations.  Using  ATS-6 
measurements  of  TEC  at  Ootacamund  in  the  Indian 
sector,  Sethia  et  al.  [1978]  report  only  one  maximum 
in  daytime  TEC  at  the  magnetic  equator. 

The  premidnight  enhancement  in  TEC  at 


±10  -15  latitude  results  from  the  postsunset  en¬ 
hancement  in  upward  drift  and  has  been  observed  at 
a  number  of  low-latitude  locations,  most  notably  at 
Ascension  Island.  Figure  13,  taken  from  the  paper  by 
Anderson  and  Klobuchar  [1983],  shows  the  large  day- 
to-day  variation  in  TEC  at  Ascension  Island  for  the 
month  of  September  1979.  Often  the  premidnight 
maximum  in  TEC  is  greater  than  the  daytime  maxi¬ 
mum  value.  Anderson  and  Klobuchar  found  that  the 
strength  of  the  meridional  neutral  wind  had  a  strong 
modulating  effect  on  the  magnitude  of  the  premid¬ 
night  enhancement.  A  strong  poleward  wind  in  the 
late  afternoon  reduces  the  TEC  enhancement  by  low¬ 
ering  the  F  layer  at  the  crests  of  the  anomaly. 

Figures  11  and  12  also  display  contours  of  TEC 
for  the  equinoctial,  solar  minimum  and  maximum 
periods  of  1975  and  1970  observed  in  the  Indian  sub¬ 
continent  region  [ Klobuchar  et  al..  1977],  The  con¬ 
tour  shapes  are  very  similar;  however,  the  observed 
magnitudes  are  higher  during  solar  minimum  and 
lower  during  solar  maximum  than  the  SLIM  calcu¬ 
lated  TEC  values.  Presumably,  this  is  because  the 
1970  and  1975  sunspot  numbers  and  F10.7  cm  flux 
values  (and  hence  ionization  production  rates)  were 
less  extreme  than  the  conditions  assumed  in  the 
SLIM  calculations. 


SUM  INDIAN  SECTOR 


LOCAL  Til*  (hr* ) 


Fig.  12.  Contours  of  calculated  TEC  values  (  x  I012  el/cnt2|, 
for  equinox,  solar  cycle  maximum  conditions,  left;  and  contours  of 
observed  TEC  values  in  the  Indian  sector  [after  Klobuchar  et  al. 
1977).  right. 
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Kig.  1 3.  Measured  values  of  TEC  versus  local  lime  at  Ascen¬ 
sion  Island  for  several  days  in  September  1979  [after  Anderson  and 
K  lohmhar,  1983] 


6  CONCLUSIONS 

A  very  brief  description  of  the  generation  of  a 
semi-empirical  low-latitude  ionospheric  model  has 
been  presented  in  which  theoretically  calculated  elec¬ 
tron  density  profiles  are  approximated  by  a  modified 
Chapman  function  given  by  (3). 

=  Nm„  exp  [<11  -  z  -  e  *')] 

where 


A 

and  Nm„  and  /im>1  represent  the  peak  electron  den¬ 
sity  and  the  altitude  of  the  F2  peak,  respectively. 
Above  the  peak  altitude,  zup  =  (h  —  hmtK)/Aup  and 
c  =  c„p  while  below  the  peak,  zl0  =  (h  -  hmtJ/Alo  and 
c  =  c,„.  The  six  coefficients  which  describe  each  pro¬ 
file  are  Nmx,  hm„,  App,  cup,  A]a  and  cl0.  Appropriate 
coefficients  which  regenerate  electron  density  profiles 
(180-1800  km),  every  half  hour  local  time  and  every 
2°  dip  latitude  between  24°N  and  24°  dip  latitude 
have  been  calculated.  Six  cases  have  been  generated, 
both  solstice  periods  and  one  equinox  period  for 
solar  cycle  minimum  and  solar  maximum  conditions. 
These  coefficients  and  resultant  electron  density  pro¬ 
files  are  provided  in  Anderson  et  al.  [1985]  and  are 
also  available  in  computer  readable  form  upon  re¬ 
quest. 


A  subroutine  which  is  able  to  quickly  generate  re¬ 
alistic  electron  density  profiles  has  a  number  of  appli¬ 
cations.  By  providing  diurnal  Ne(h)  and  TEC  values 
at  any  specified  dip  latitude  or  slant  path,  SLIM 
could  be  used  to  examine  propagation  effects  deter¬ 
mined  by  ionospheric  structure.  SLIM  has  several 
possible  applications  related  to  so-called  “active  ex¬ 
periments"  where  ambient  ionospheric  conditions 
need  to  be  specified  for  preexperiment  planning  and 
postevent  analysis.  In  addition,  such  a  model  would 
yield  flux-tube  integrated  Pedersen  conductivity  and 
electron  content  values  which  are  important  in  deter¬ 
mining  and  predicting  low-latitude  instability  growth 
rates.  Finally,  the  program  would  be  capable  of 
supplying  airglow  intensities  both  in  the  vertical  and 
at  any  slant  viewing  angle  desired. 

One  final  note  of  caution,  since  the  theoretical  cal¬ 
culations  involve  only  O*  ions,  SLIM  is  basically  an 
F  region  ionospheric  model  and  does  not  give  realis¬ 
tic  electron  densities  below  the  F  layer  where  NO* 
and  02  may  dominate  or  the  topside  where  H  +  may 
be  the  dominant  ion.  Daytime  and  nightime  electron 
densities  in  the  topside  region  are  typically  around 
104  el/cm3  with  very  laTge  scale  heights  reflecting  the 
low  mass  of  the  H  *  ion.  During  the  day,  bottomside 
electron  and  ion  densities  in  the  110  to  180  km  alti¬ 
tude  range  are  1  to  3  x  10s  el/cm3  and  depend  on 
solar  zenith  angle.  At  night,  the  low-latitude  bottom- 
side  NO*  density  profiles  are  highly  variable  with  F, 
region  densities  varying  between  1  x  105  and  2 
x  104  el/cm3.  It  may  be  necessary,  therefore,  for  the 
users  to  supply  their  own  electron  density  specifi¬ 
cation  in  the  ionospheric  regions  where  SLIM  is  not 
valid. 


APPENDIX:  ANALYTIC  WIND  EQUATIONS 
Equinox 


V0  =  8.61  x  I0'3»fF  — 


ft 

io\ 


2  n(  |0|\ 

V2A  =  20  co*  —  ( f  —  18  -  — ) 

24  V  7  / 

(  n  10  \  In  /  |0|\ 

V,  2  =  50  sin  ( - »  —  ft  )  cos  —  ( (  -  9 - ) 

\180  3  /  12  V.  12/ 

2  n 

K,  =  0.380  cos  —  (f  -  8.2  +  0.079 1 0|) 

8 

I'tot  =  K}  +  -I  U,  2  +  Fg 
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The  6  is  latitude  in  degrees. 

December  solstice 

V„  =  3.116  x  10 -**e\0  +  45)  -  49 


2n  ( 

^*  =  20„cos  —  If  -  18 - j-\  0J4  =  0+1O 


V\  2  =  50  sin 


fJL.i® 

\180  3 


2n 


V,  =  0.380  cos  —  (/  -  8.2  +  0.079 10|) 

8 

^TOT  -V0+Vl*+V,2+ 

June  solstice 

V0  =  -3.116  x  10  4  «  0*(45  -  0)  +  49 
2n  /  |024|\ 

U24  =  2024  cos  -  /  /  -  18 - O2.  =  0- 

(  n  10  \  2n  (  |0|\ 

y,  2  =  50  sin  ( - «  —  0 )  cos  —  |  r  —  9 - ) 

\  180  3  /  12  \  12/ 


yH  =  0.380  cos  —  (I  -  8.2  +  0.079|0|) 
8 

where 


10 


I'TOT  —  4-  ^ 2 *  7"  1^12  4-  (m/S) 

and  t  is  hours,  local  time;  9  is  geographic  latitude, 
positive  northern  hemisphere. 
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